arXiv:hep-ph/9701235vl 7 Jan 1997 


SNUTP 96-094 


Color-octet mechanism in the inclusive D—wave 
charmonium productions in B decays 


Pyungwon Ko^’* 

^Department of Physies, Hong-Ik University, Seoul 121-791, Korea 


Jungil Lee^’'^ and H.S. Song^d 
‘^Department of physics and Center for Theoretical Physics 
Seoul National University, Seoul 151-742, Korea 
(February 1, 2008) 


Abstract 

The inclusive L)—wave charmonium production rates in B decays are consid¬ 
ered in the Bodwin-Braaten-Lepage (BBL) approach. We find that the color- 
octet subprocesses B —> cc(^5|®Ar -|- s,d, followed by the transition 

cc{^s[^\ or ^Pj^'^) —>■ ^Dj, strongly dominate over any other subprocess, due 
to the large Wilson coefficient for the AB = 1 effective lagrangian. Assum¬ 
ing that the numerical values of the matrix elements (0|O^^‘^(^S'j*^)|0) and 
(0|C>^^"'(^pj^^)|0) are the same order of magnitudes with the (0|C>^ (^5®)|0), 
we find that the ^’^D 2 can be observed at future B-factories. 
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1 . Estimates of S— and P—wave charmonium production rates in B decays were one of 
the earliest applications of the heavy quarkonium physics in the framework of perturbative 
QCD (PQCD) [|^]. In the case of P—wave charmonia, it had long been known that their decay 
and production rates in the color-singlet model have infrared singularities, which signaled 
the failure of the factorization procedure in terms of a single wavefunction of the cc pair in 
the color-singlet P—wave state. Recent development in Nonrelativistic QCD (NRQCD) as 
an effective field theory of QCD helped such inconsistency be resolved including higher 
Fock states of the P—wave charmonia in a systematic manner in NRQCD [^]. 

On the contrary, the S'—wave charmonium productions in B decays have no such infrared 
divergence problems. However, the prediction based on the color-singlet model falls short 
of the data by a factor of ~ 3. This can be partly overcome by including the color-octet 
mechanism, since the relevant Wilson coefficient is a factor of ~ 30 larger than that of the 
color-singlet case p. Also the polarization of J/V’ produced in B decays depend sensitively 
on the color-octet contributions 0. 

In this work, we extend previous works i § i to the case of the P—wave charmonium 
productions in B decays. In various potential models for the charmonium family, P—wave 
charmonium states are predicted in the mass range of 3.81 — 3.84 GeV []^. Although the 
masses are above the PP threshold, they are below the DD* threshold. Also, two states 
= 2—) and ^P 2 (J^^ = 2-+) are forbidden to decay into PP because of parity 
conservation. So these two states are predicted to have narrow decay widths, and we can 
easily tag these two states through dominant decay channels (to be discussed below). There¬ 
fore, P—wave charmonium states (at least for J = 2) is certainly accessible at P-factories 
where one expects to have ~ 10® B decays per year. Thus B factories can serve as a place 
for studying new charmonium states below/above the PP threshold, including some hypo¬ 
thetical states like hybrids or Xcj{‘^P) state |p, in addition to its main role for measurements 
of CP asymmetries in B decays. 

2. The AB = 1 effective Hamiltonian for b ccq ( with q = d,s) is written as [Q 
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( 1 ) 


He,, = ^ V^v:, 


2C+ - C- 


c7m( 1 -75)c g7^(l - 75)6 


+ (C+ + C_) C7m(1 - 75)T“c qY{l - -f,)T%] , 


where are the Wilson coefficients at the scale /r ~ M^. We have neglected penguin 
operators, since their Wilson coefficients are small and thus they are irrelevant to our case. 
To the leading order in as{Mb) and to all orders in as{Mb) ln{Mw/Mb), the above Wilson 
coefficients are [Q 


C+{Mb) ^ 0.87, C_{Mb) ^ 1.34. (2) 

At the parton level, the B decay into charmonium state occurs through 

b + s,d^ ^ ^ 3 ^ 

Let us hrst estimate the velocity scaling of the various subprocesses. From the effective 
Hamiltonian given in Eq. O. one can get the amplitude for b —> + s, d by 

the Taylor expansion of the amplitude with respect to the relative velocity (n) of the heavy 
quark inside the heavy quarkonium and projection over the specihc partial wave state of 
(cc). Expanding the physical state in terms of higher Fock states as following, 

\^Dj) = 0{l)\ccCD^J^)) + 0{v^)\cc{^PP)g)+0{v^)\^^^^ 

+ 0{v^)\ccCDf)gg) + 0{v^)\cc{^s[^'^^)gg) + 0{v^)\ccCD^2^)g) + ..., (4) 

= 0{1)\cc{^D^2^)) + 0{v^)\ccQp[^^)g) + O(n^) |cc(^Ff 

+ 0{v‘^)\ccQDf^)gg) + 0{v‘^)\cc{^S'i'^^)gg) + 0{v‘^)\cc{^Df)g) + ..., (5) 

we can End the velocity scaling of the soft process. Note that the chromoelectric El and the 
chromomagnetic Ml transitions accompany with the extra factors of v and respectively, 
in the amplitude for the soft transition, —»• In Table |, we present 

the velocity scalings of the amplitude in NRQCD for various channels. According to Table |, 
the following processes are leading order in velocity scaling : 
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+ 

s, d 
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s, d 

cc) 
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^ ^Dj + X, 

^ ^D2 + X. (7) 

Note that in the hard process amplitude for B —> j + s, d, the spectroscopic 

states do not contribute in the leading order. Since the Wilson coefficient for the 

color-octet subprocess is about 30 times larger than that of the color-singlet subprocess in 
the (amplitude)^ level P], we hnd that the following subprocesses (which are 0{v^) in the 
amplitude level) dominate over all the others : 


b —>• {cc)^Pj^^ + s,d 
b —>• (cc)^S'l®^ -|- s, d 

b —>■ (cc)^S'® + s,d 


^Dj + X, 
^D2 + X, 


( 8 ) 

(9) 


Using the results in Ref. we immediately obtain the inclusive charmonium pro¬ 

duction rate in B decay as 


T{B 


^Dj + X) ~ V{B ^ {ccfsf'^ + s,d^ ^Dj + X) 
+ T{B ^ {ccfB^ + s,d^ ^Dj + X) 


= {C+ + C.f 1 + 


8 M 2 

'W 


X 




r(p 


2 M 2 

lP 2 + X) ~ r(P ^ + S, d ^ ^P 2 + X) 

3(0|0'^2(iS'f^)|0) 


= (U+ + C'_)^rc 


where 


fo = IWfcP 




2 M 2 


( 10 ) 


( 11 ) 
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( 12 ) 
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In order to evaluate the branching ratios numerically, we need to know three nonperturbative 
matrix elements, |0), (0|0^^2(^S'o^^)|0). 

In connection with this, we note that Qiao et. al. have recently shown that the color- 
octet production processes —>• ^Dj{cc)qq have distinctively large branching ratios and it 
is of the same order as that of the J/'ifj production in decays, assuming 

5 

= X 4.6 X 10-3 GeV3, for J = 1, 2, 3. (13) 

5 

Similarly, we can estimate the numerical value of the matrix element (0|(!1^'^‘^(3 Pq 3^)|0) as 

5 

= X 1.3 X 10-2 GeV^ for J = 1, 2, 3, (14) 

0 


using the numerical value which was maximally allowed in Ref. []^. This value is about an 
order of magnitude larger than that deduced from the velocity scaling. For the case of ^P 2 , 
we use the numerical value of by imposing the approximate heavy quark 

spin symmetry as 




( 16 ) 


With the numerical values as 


Mfe = 5.3 GeV, = 1.5 GeV, (16) 

a,(2MQ = 0.253, (17) 

(0|O'^2(3Ff^)|0) = 4.6 X 10-3 GeV3, (18) 


we estimate the branching ratios as 

. 2J+ 1 \ 5.2% when ( 0 |O'^R 3 p( 8 ) 6 |Qi _ ^ 3 ^ lo-SQeV^ 

P(P^ 3p^p^)^^T I ojl / ^ 

^ [ 1 . 1 % when ( 0 |O ^ 2 ( 3 P(S 3 ))|Q^ _ 3 3 ^ lo-SQeV^ 

P(P ^ ^D2 + X) = 0.67%. (20) 


5 








Before closing this subsection, we have to mention on theoretical uncertainties in our 
estimates for B{B —^D + X). The hrst comes from the assumptions, (13)-(15), which 
originate from our ignorance of nonperturbative physics in QCD. Up to now, there are no 
known processes which depend on these three parameters. So we relied upon the naive 
guess based on the velocity scaling rule and heavy quark spin symmetry. Another intrinsic 
uncertainty comes from the approach adopted here and in Ref. [§]. Namely, the parton 
picture adopted here for the case of D—wave charmonium productions in B decays may not 
be so good as in the case of the S'—wave charmonium productions in Ref. p because the 
phase space available in the former is much less than the latter. There are also uncertainties 
coming from the quark masses Mf,, and the relation between and the mass of the 
D—wave charmonium state. In order to minimize these last uncertainties, we have used the 
prescriptions made in Ref. [Q. Despite of all of these uncertainties, we still anticipate our 
results to be correct estimates of an order of magnitude for the branching ratios for the 
decays B ^ D + X. 

3. In order to identify the desired D—wave charmonium state produced in B decays, one 
has to know its decay channels. As mentioned earlier, J = 2 states (^'^■'■^^ 2 ) are predicted to 
have narrow decay widths of order of 300 — 400 keV. So these states are easier to reconstruct 
in the experiments, and hereafter we consider only J = 2 D—wave states, and try to find 
out which mode is best to identify the J = 2 D—wscve charmonium state in B decays. 

First consider the spin-triplet D—wave state (= 2 ). The main decay modes of this 
state are the El transitions into Xcj^j 2 (lP) states, the hadronic transition into J/tf + tttt, 
and the decay into light hadrons via cc annihilation into three gluons. The decay widths 
for each process is rather model dependent, but it is enough to quote some nominal value 
for each mode, in order to have some idea about the branching ratio for each channel. We 
adopt the numbers given in Ref. : 

5('^2^Xci + 7) =0.64, 

5('/^2^Xc2 + 7) = 0.15, (21) 
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B{^D2 —> J/ijj + vr’''7r ) = 0.12. 

Combining the branching ratios for Xcj^i 2 ^^^) J1'^ + 7 J/'^ to 

B{^D 2 J/iJ + vr+TT-) X Bi^J/ij /r+/i-) = 7 x 10-^ 

B{^D 2 Xci + 7 ) X B{xci W + 7 ) X BiJ/i) /i+/i”) = 10.4 X 10“^, (22) 

BCD 2 Xc2 + 7 ) X B{xc2 W + 7 ) X BiJ/i) = 1.2 X 10"^. 


The last decay channel through the intermediate 7^2 (1-P) state is negligible. For the radiative 
decay channels, one needs a good photon detector. Otherwise the hrst option through 
^D 2 J/ij + TT+TT” would be the best to reconstruct ^D 2 state at B factories. Combining 
the branching ratios for the cascade decays, we get 


B{B ^ ^D2 ^ /i+/i-) + TT+TT-) 

[ 3.7 X 10-^ when (0|0'^2(3Po(®^)|0) = 1.3 x 10-2GeV^ 

~ < (23) 

I 7.5 X 10-5 when (0|O"^2(5P({^^)|0) = 1.3 x lO-^GeV®. 

Gomparing these numbers with those for B —>■ J/'0(or ijj') + X, 


B{B 

J/ij + x- 

+X) =4.7x 10-^ 

(24) 

B{B 

^ + X - 

+X) = 2.0 X 10-^ 

(25) 


we hnd that the P—wave charmonium productions in B decays are not so rare, and can 
be within the reach of B-factories. There would be several tens to a few hundreds’ ^P 2 ’s 
in 10® B meson decays into B —> ^D 2 —>■ J/'ip'K^TT~ with J/-^ —> depending on the 


numerical values of (0|O^^^(5S'|®^)|0) and (0|0^^^(^Po®)|0). 
Next consider the spin-singlet D—wave state = 2 '“). 


for this state is 11 


One possible decay channel 


^he('Pi)+7) = 0.80 

B{KCPi) J/ij + 7r°) = 0.005 (26) 

B{J/'ip —> iX) = 0.06. 
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Here one has to tag three photons, one from the first chain, and the other two from the 7r° 
decay. The product of three branching ratios in this channel is 2.4 x 10“^, so that we expect 

B{B ^ ^ 7 + 7r° + ~ 1-6 x 10“®, (27) 

corresponding to ~ 160 such cascade events in 10®i? decays. 

4. In conclusion, we have considered the wave charmonium productions in B decays. 
The dominant contributions come from the color-octet (cc)^S'® and states for ^Dj 

and (cc)^S'® state for ^D 2 production, respectively, because of the large Wilson coefficient, 
(C+ -I- C_) in Eq. (|l]). Assuming = {0\O^ x (2J -|- l)/5 and 

(0|0'^^(3Po^^^)|0) = (0|0’^'(3Po^®^)|0) X (2J + l)/5, we obtain B{B ^D 2 + X) ^ 1.1% - 
5.1%. Assuming B{^D 2 —>■ J/ip + ^ 10% and using B{J/ip —> = 0.06, one can 

hnd 75 to 370 ^D 2 ^s in the 10® B decays into B —>• ^D 2 + X followed by ®P 2 J/ip + vr+vr” 
and J/tp ^ Pl~ with I = /x. This is quite accessible in P-factories where 10® B meson decays 
will be accumulated per year [[ Since the P—wave charmonium states are not easy to hnd at 
other colliders, it is recommended to search for these states as well as others such as Xcj (2P) 
and hybrid charmonia from B decays. Finding such states at the level predicted in this work 
will be another test of the factorization based on the NRQCD, including the color-octet 
mechanism for the heavy quarkonium productions. Also it will deepen our understanding of 
the heavy quarkonium spectroscopy based on the color-singlet pair and exotic hybrid states. 
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TABLES 



subprocess 


L 

E 

M 

total (^-^-+-^+2^) 

b- 

^ {cc)^+ s, d ^ 

^Dj + X 

2 

0 

0 


b- 

> {cc)^+ s, d ^ 

^Dj + X 

2 

2 

0 


b- 

> (cc)^ + s, d ^ 

^Dj + X 

2 

0 

1 


b^ 

-|- s, d — 

> ^Dj + X 

1 

1 

1 

v'^ 

b- 

■> {cc)^pf'^ + s,d ^ 

^Dj + X 

1 

1 

0 


b^ 

■ (cc)^5q^’®^- 1-s,d ^ 

► ^Dj + X 

0 

2 

1 


b^ 

• -|-s, d ^ 

► ^Dj + X 

0 

2 

0 


b- 

(cc)^L>2^^ -|- s, d ^ 

^D2+X 

2 

0 

0 

v‘^ 

b- 

(cc)^L>2*^ -|- s, d ^ 

^D2+X 

2 

2 

0 


b- 

{cc)^D^P + s,d ^ 

^D2+X 

2 

0 

1 


b^ 

• {cc)^Pj^’^^ + s,d — 

> ^D2 + X 

1 

1 

1 


6 ^ 

► {cc)^Sq^’^^ + s,d — 

> ^D2+X 

0 

2 

0 

v‘^ 

6 ^ 

► -|-s, d — 

> ^D2+X 

0 

2 

1 

v'^ 


TABLE I. The velocity scalings of the amplitudes for the the various subprocesses relevant 
to B ^ ^Dj + X and B ^D 2 + X. L is the number of derivatives in the NRQCD 4-quark 
operators {O ^and the E and M are the numbers of the chromoelectric El and the 
chromomagnetic Ml transitions in order to reach the color-singlet physical H—wave state. 
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